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Abstract 

A Bianchi type-I cosmological model in the presence of a magnetic 
flux along a cosmic string is investigated. A nonlinear spinor field is 
used to simulate the cosmological cloud of strings. It is shown that 
the spinor field simulation offer the possibility to solve the system of 
Einstein's equation without any additional assumptions. It is shown 
that the present model is nonsingular at the end of the evolution and 
does not allow the anisotropic Universe to turn into an isotropic one. 
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1 Introduction 

Though the present day Universe is well described by an isotropic and homo- 
geneous Friedmann-Robertson- Walker (FRW) model, there are serious theo- 
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retical arguments about the existence of an anisotropic phase in the evolution 
of the Universe pQ. These ideas were further supported by the observational 
data from COBE (Cosmic Background Explorer) and WMAP (The Wilkin- 
son Microwave Anisotropy Probe) where small anisotropy in the microwave 
background radiation was found. These lead many cosmologists to consider 
the anisotropic model to describe the initial phase of the evolution which 
eventually decays into an isotropic FRW one. 

After the famous paper by A. Guth [2], scalar field was hugely used in 
simulating different cosmological models. But the question occurs, if other 
fields can contribute to the evolution of the Universe. As an alternative the 
spinor field was used due to its sensitivity to the gravitational one [21 IH 
El El El El El ID]- The nonlinear spinor field proved to be able to describe 
different cosmological models. For example it was shown that a suitable 
choice of nonlinearity (i) accelerates the isotropization process, (ii) gives rise 
to a singularity-free Universe and (iii) generates late time acceleration. 

Moreover, it is shown that a nonlinear spinor field can be used to simulate 
a perfect fluid from ekpyrotic matter to phantom matter [TTl [T21 fT3] . 

At the same time the string cosmological models have been used in at- 
tempts to describe the early Universe and to investigate anisotropic dark 
energy component including a coupling between dark energy and a perfect 
fluid (dark matter) [HI [15]. Cosmic strings are one dimensional topological 
defects associated with spontaneous symmetry breaking in gauge theories. 
Their presence in the early Universe can be justified in the frame of grand 
unified theories (GUT). 

On the other hand, the magnetic field has an important role at the cosmo- 
logical scale and is present in galactic and intergalactic spaces. Any theoret- 
ical study of cosmological models which contain a magnetic field must take 
into account that the corresponding Universes are necessarily anisotropic. 
Among the anisotropic spacetimes, Bianchi type-I space (BI) seems to be 
the most convenient for testing different cosmological models. 

The object of this paper is to investigate a BI string cosmological model in 
the presence of a magnetic flux. For this purpose we use the nonlinear spinor 
field simulation as was described recently in a number of papers [HI [T^l IT3] . 
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2 Basic equations 



We study the evolution of the Universe in presence of a cosmic string and a 
magnetic flux in the framework of a BI anisotropic cosmological model. For 
a BI spacetime the line element is given by 

ds^ = {dtf - al{t){dx^f - al{t){dx^f - al{t){dxY . (1) 

There are three scale factors (i = 1, 2, 3) which are functions of time t only 
and consequently three expansion rates. In principle all these scale factors 
could be different and it is useful to express the mean expansion rate in terms 
of the average Hubble rate: 

H=^-(^ + ^ + ^), (2) 



3 Voi 02 cts 

where over-dot means differentiation with respect to t. 

Einstein's gravitational field equations, corresponding to the metric ([1]) 
have the form 

02^03^ 02 03 ^ ^^1^ ^3^^ 

02 O3 O2 O3 

03 Oi 0-3 di 2 /QU\ 

\ \ = , (3b) 

O3 Oi O3 Oi 

^ + ^ + ^^ = KTi, (3c) 

Oi O2 Oi O2 



Ol O2 O2 O3 O3 Ol Q 

\ \ = kTq , (3d) 

Ol O2 O2 O3 O3 Ol 

where k is the gravitational constant. The energy momentum tensor for a 
system of cosmic string and magnetic field in a comoving coordinate is given 
by 

T; = puX - Aa;^a;^ + E; , (4) 

where p is the rest energy density of strings with massive particles attached 
to them and can be expressed as p = Pp + A, where Pp is the rest energy 
density of the particles attached to the strings and A is the tension density of 
the system of strings |T6l ITTl |18] which may be positive or negative. Here Ui 
is the four velocity and Xi is the direction of the string, obeying the relations 

UiU^ = —XiX^ = 1, Uix"^ = . (5) 
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In dlj E^i, is the electromagnetic field [TP] and in what follows we shall 
choose the string and also the magnetic field along direction. In our model 
the electromagnetic field tensor F°'^ has only one non- vanishing component, 
namely 

F23 = h, (6) 

where h is presumed to be constant. For the electromagnetic field one 
gets the following non-trivial components 

where /i is a constant characteristic of the medium and called the magnetic 
permeability. Typically p, differs from unity only by a few parts in 10~^ 
(/i > 1 for paramagnetic substances and /i < 1 for diamagnetic) . To simplify 
the notation, we include in a constant /3 the value of the electromagnetic 
field, h, and the magnetic permeability, /2. 

Using comoving coordinates we have the following components of energy 
momentum tensor [2U] : 

TS-P = TI~X = -T| = -T| = (8) 
where we introduce the volume scale of the BI space-time [1] 

r = = 010203 , (9) 

which is connected with the Hubble rate ([2]), namely ^ = 3H . 
In view of T2 = from ([8]) one finds 

02 = 03/^exp(^X j , (10) 

with D and X some integration constants. 

On the other hand, summation of Einstein equations (15a|l . ( l3b|) . (|5c]) and 
three times ( l3d|) gives: 

£ = iK(3p + A + ^=^). (11) 

Taking into account the conservation of the energy-momentum tensor, 
i.e., T^.^ = 0, after a little manipulation of ([8]) one obtains [2T| [22]: 

7" O1 , . 

P + -P- —X = 0. 12 

r Oi 
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It is customary to assume a relation between p and A in accordance with 
the state equations for strings. The simplest one is a proportionality relation 

p = aX . (13) 
Among the most usual choices of the constant a we mention the following: 

{1, geometric string 

1 + uj > 0, p string or Takabayasi string (14) 

— 1, Reddy string 

It is possible to consider a more general barotropic relation 

p = p(A) . (15) 

than the linear relation f|T3l) . subject to the restrictions imposed by the energy 
conditions. The weak energy condition as well as the strong one require p > X 
with A > or p > with A < and the dominant energy condition implies 
p > and p2 > [l6j. 



3 Spinor field approach 

Recently it was shown that a nonlinear spinor can be used to simulate differ- 
ent types of perfect fluids including those called ekpyrotic, phantom matter 
and dark energy [HI [121 113] • Here we show that it is possible to describe 
cosmic strings in terms of spinor fields as well. 



3.1 Spinor field simulation 

We shall simulate the cloud formed by massive cosmic strings with particles 
attached along their extensions with a nonlinear spinor field described by the 
Lagrangian: 



— miptp + F , 



(16) 



with F being some arbitrary function of the scalar S = '4)il). For the simula- 
tion of the present cosmological model in which the anisotropic scale factors 
Oj are functions solely of time it is adequate to assume that the spinor field 
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depends on t only. Then the corresponding components of energy-momentum 
tensor take the form 

T° = mS-F, (17a) 
r/F 

V = S^-F ,z = l,2,3. (17b) 

In what follows we describe the energy density p of the string by Tq and 
the tension density A by in agreement with ([8]). Inserting fll7ap and (117bp 
into (fTSll we find 

S^-il--)F-^S = 0, (18) 
ab a a 

with the solution 

F = + mS . (19) 

Here Ci is an integration constant. The positivity of Tq imposes some re- 
striction on Ci, namely Ci < 0. Setting Ci = —u we find the cosmic string 
can be described by the spinor field Lagrangian 



^sp-- 



- z/S^"-^)/" , (20) 



remarking the disappearance of the mass term [TT 
With these preparatives we get 



p = T° = z/5("-^)/", (21a) 
X = T^ = (21b) 



Variation on (U61) with respect to ip and ip gives spinor field equations with 
nonlinear terms. On the other hand from the spinor field equations for 5* one 
finds 

S + -S = 0, (22) 

T 



with the solution 



S=^, (23) 
r 



Co being a constant. 

Taking into account this simple behavior of S we have finally 



a-l 



p = Rt~— , (24a) 

R a-l , , , 

A = -r"— , 24b 
a 



with the constant R = z/Cq " . 

Using the above formulas for p and A, from f lT2|) we can determine the 
anisotropic factor ai: 

ai = A,T, (25) 

Ai being a constant of integration. On the other hand, from equations ([9]) 
and ffTOj) we obtain 

ro" fx f dt'\ , , 

and 

1 / X fdt'\ 

3.2 Asymptotic behavior 

In what follows, we study the equation for r in details. Using the above 
expressions for p, A, ai we get from flTT]) : 

r = Rt^ + — — . 28) 

2 a 2 ^ ^ 

We can evaluate the derivative of r with respect to t which leads finally to a 
solution in quadrature 

dr 

t + to, (29) 



^/[nR{3a + l)/(a + l)]r("+i)/" + [Kf3^Al/2]T^ + C 

C and to being some integration constants. 

In spite of the fact that this equation cannot be explicitly solved, the 
asymptotic behavior of the solutions for t — )• oo could be found. 

3.2.1 Case I 

In most cases, fl2^ provides a standing expansion of the volume scale of the 
BI spacetime for t growing. 
Indeed, for 

a + 1 , ^ 

—— < 2 , 30 

a 
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i.e. a > 1 or a < 0, the term with at the denominator of the 1. h. s. of 
f l29p is dominant and we get an exponential behavior 

T oc exp t , (31) 

for large t. In this case, only the anisotropic factor ai presents an exponential 
increase for t — )■ oo, while the scale factors 02 and 03 tend to constants. 

Concerning the asymptotic behavior of p and A we infer from fl24ap , fl24bp 
that they tend to zero as 

P^^°^ 

exp [—t) 



3.2.2 Case II 

For < a < 1, there is no solution with a r — 00 behavior for t — 00. In 
this case from (129|) one finds: 

r 2q oc t. 

Since a — 1 < 0, r cannot tend to infinity as t — )■ 00. Therefore the model 
does not admit a consistent solution for < a < 1 in agreement with the 
discussion of the general barotropic equation f lTSj) from Section [2j 

Taking into account the absence of consistent solutions for < a < 1 
as it was shown above, for the present spinor simulation, in what follows we 
shall refer only to the situations described in subsection 13.2.11 (Case I). It is 
worth also mentioning that the present model does not support in any case 
a vanishing of r for t — )■ 00. 



3.3 Numerical simulations 

In this subsection we graphically illustrate the evolution of energy density p, 
volume scale r and metric functions ai, a2 and 03. In doing so we rewrite 
the equation for r f l28|) in terms of the Hubble parameter H and introduce a 
new function T: 

f = 3Ht, (33a) 

E = ~3H^ + - —Rt— + ---^, (33b) 

b a D 

T 

T = -. (33c) 
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We also rewrite the metric functions in terms of T, which now read 

a^^A.r, a.-^T'"^ a,^ ^r'^'K (34) 

Let us now numerically solve the system fl33|) . In doing so for simplicity 
we set K = 1, = 1, D = 1 and X = 1. In Fig. [T]we plot the evolution of 
energy density p for different values of a, namely for a = —1 and a = 1.5. 
Evolution of r corresponding to these values of parameters in shown in Fig. 
121 For simplicity we also set the initial values of r, H and T to be unity. 




t 



Figure 1: View of energy density Figure 2: Evolution of the Uni- 

of the cosmic string for different verse corresponding to the energy 

value of a. densities given in Fig. [1] 

The numerical simulations support the behavior described in (1311) . (132|1 . 
From Fig. |2] one finds at first, in the case of a negative a the Universe 
expands slower than it does as for a positive a, though in both cases we have 
an exponential growth for large time. 

In Figs. [3] and H] we plot the behavior of metric functions for positive 
and negative values of a. As is seen from the figures, with the expansion of 
the Universe ai increases exponentially, while 02,03 tend to constants. Thus 
we see that introduction of cosmic string does not allow isotropization of 
initially anisotropic space-time. In the following subsections we discuss the 
singularity problem and isotropization process in detail. 
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Figure 3: Evolution of metric func- 
tions for a = 1.5. 



Figure 4: Evolution of metric func- 
tions for a = — 1. 



3.4 Singularities 

The next task is to investigate the space-time singularities in the present 
model. Let us rewrite the metric functions in the following form: 



tti = CiT^' exp ( Yi 



dt' 



(35) 



where Ci = Ai, C2 = /d/A, C3 = l/^DM, = +1, N2 = N3 = 0, Y, = 
0, Y2 = X/2, and Y3 = —X/2. Then the first and the second derivatives of 
the metric functions take the form: 



Oh 
a,- 



T T 



a,- r Vr/ t'^ 



(36a) 
(36b) 



We study the singularities analyzing the regularity properties of the Krets- 
chmann scalar which for the metric ([1]) reads 



X 



+ 



+ 



ai a2 
ai a2 



a2 as 
02 cts 



03 oi- 
03 fli - 



• (37) 
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Evidently, a singularity can occur when some or all scale factors tend to 
zero or infinity. We shall follow the criteria given in [23] which states: 

(i) t finite, some Oj — t- oo. // at least one scale factor becomes infinite at 
finite t, it is a curvature singularity. 

(ii) t finite, some Oj — i- 0. // more than one scale factor turns to zero at 
finite t, it is a singularity. If only one scale factor is zero at finite t, the 
spacetime can be nonsingular. 

(iii) t — )■ oo, some a, — )• oo. Such an asymptotic can only be singular if at 
least one scale factor grows faster than exponentially: 

aiit) ^ exp k = const. > 0. 

(iv) t — )■ oo, some Oj — )■ 0. Such an asymptotic can only be singular if at 
least one scale factor vanishes faster than exponentially: 

aj(t) < 0(exp A; = const. > 0. 

Looking at the asymptotic behaviors described above (Section 13. 2p we 
conclude that in the spinor field approach of the cosmic strings the present 
model is nonsingular at the end of the evolution. No scale factor of the BI 
metric presents a growth faster than exponentially or vanishing faster than 
exponentially for t — > oo. 



3.5 Isotropization 

Since the present-day Universe is surprisingly isotropic, it is important to 
see whether our anisotropic BI model evolves into an isotropic FRW model. 
Isotropization means that at large physical times t, when the volume factor 
r tends to infinity, the three scale factors ai{t) grow at the same rate. Two 
wide-spread criteria of isotropization read 

i=l 

T? = ^AH^^O. (38b) 

Here A and are the average anisotropy and shear, respectively, while 
Hi = di/ui is the directional Hubble parameter evaluated in ( I36ap . We also 
investigate the isotropization condition proposed in 

a 



— )■ const , (39) 

t— s-oo 
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where a{t) = the average scale factor. Provided that condition fl39p 

is fulfilled, by rescaling some of the coordinates, we can make Oj/a — )■ 1, and 
the metric will become manifestly isotropic at large t. 

In order to study the isotropization process in the present model, we 
evaluate from ( 135|) 

Taking into account the value of Ni and Yi we see that Qi/a do not tend to 
constants as t — ?■ oo. For example, in the generic case with the exponential 
expansion of the volume of the Universe (13T|) . ai/a — )■ oo, 02/0 — )■ and 
fls/a — i- 0. So in the case considered, no isotropization process takes place. 
Indeed, at the early stage of evolution, where r — )■ we assume Oi — )• 00, 
a2 — )• and 03 — )■ 0, that is at this stage the Universe looks like a one- 
dimensional string. In the asymptotic region where t ^ 00 and r — 00 we 
have ai — 00, while a2 and 03 evaluating to finite values. 

Let us also examine the possibility of an isotropization process using 
conditions ( 155]) . In view of fl5^ from f l55]) in the generic case f l5T]) one finds: 

A = 2, (41a) 
= i. (41b) 

As one sees, both A and are some positive defined, non vanishing quan- 
tities. Thus, in general, the present model does not undergo an isotropization 
process, the anisotropic feature of the Universe from the early stages being 
retained. 

4 Conclusions 

We have studied the evolution of a anisotropic universe given by a BI cos- 
mological model in presence of cosmic strings and magnetic field. In doing 
so we exploit the spinor simulation of cosmic strings that allows us to solve 
the system of Einstein's equations. Exact solutions have been supplemented 
with some numerical evaluations. It is shown that the presence of cosmic 
strings does not allow the anisotropic Universe to evolve into an isotropic 
one. 



12 



Acknowledgments 

The authors gratefully acknowledge the support from the joint Romanian- 
LIT, JINR, Dubna Research Project, theme no. 05-6-1060-2005/2010. M.V. 
is partially supported by the GNOSIS Program IDEI - 571/2008, Romania. 

References 

[1] C.W. Misner, The Astrophys. J. 151, 431 (1968). 
[2] A. Guth, Phys. Rev. D 23, 347 (1981). 
[3] M. Henneaux, Phys. Rev. D 21, 857 (1980). 
[4] B. Saha, Phys. Rev. 64, 123501 (2001). 

[5] G. Armendariz- Picon and P.B. Greene, General Relat. Grav. 35, 1637 
(2003). 

[6] B. Saha and T. Boyadjiev, Phys. Rev. D 69, 124010 (2004). 

[7] V.O. Ribas, F.P. Devecchi, and G.M. Kremer, Phys. Rev. D 72, 123502 
(2005). 

[8] B. Saha, Physics of Particles and Nuclei 37. Suppl. 1, S13 (2006). 
[9] B. Saha, Phys. Rev. D 74, 124030 (2006). 
[10] B. Saha, Phys. Rev. D 69, 124006 (2004). 

[11] V.G. Krcchet, M.L. Fel'chenkov and G.N. Shikin, Grav. Cosmology 14, 
292 (2008). 

[12] B. Saha, Spinor model of a perfect fluid , arXiv: 0901.1387 [gr-qc]. 

[13] B. Saha, Spinor model of a perfect fluid: examples, arXiv: 0902.2097 
[gr-qc]. 

[14] T. Koivisto and D.F. Mota, Phys. Lett. B 644, 104 (2007). 

[15] T. Koivisto and D.F. Mota, Phys. Rev. D 75, 023518 (2007). 
[16] P.S. Letelier, Phys. Rev. D 28, 2414 (1983). 

13 



[17] A. Pradhan, A. K. Yadav, R. P. Singh and V. K. Singh, Astrophys. Space 
Sci. 312, 145 (2007). 

[18] G.S. Khadekar and S.D. Tade, Astrophys. Space Sci. 310, 47 (2007). 

[19] A. Lichncrowicz, Relativistic Hydrodynamics and Magnetohydrodynam- 
ics, (Benjamin, New York, 1967). 

[20] B. Saha, Astrophys. Space Sci. 299, 149 (2005). 

[21] B. Saha and M. Visinescu, Astrophys. Space Sci. 315, 99 (2008). 

[22] B. Saha, V. Rikhvitsky and M. Visinescu, Cent. Eur. J. Phys. 8, 113 
(2010). 

[23] K.A. Bronnikov, E.N. Chudaeva and G.N. Shikin, Class. Quantum Grav. 
21, 3389 (2004). 



14 



